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Abstract 

We present a construction of fermionic operators in 3+1 dimensions in terms 
of bosonic fields in the framework of QED4. The basic bosonic variables are the 
electric fields Ei and their conjugate momenta Ai. Our construction generalizes the 
analogous constuction of fermionic operators in 2+1 dimensions. Loosely speaking, 
a fermionic operator is represented as a product of an operator that creates a 
pointlike charge and an operator that creates an infinitesimal t'Hooft loop of half 
integer strength. We also show how the axial U{1) transformations are realized in 
this construction. 
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The problem of bosonization in dimension higher than two is a long standing one. 
Many attempts have been made to generalize the 1+1 dimensional Mandelstam [jl| con- 
struction to higher dimensions 0. However, they all have an unattractive feature that 
the bosonized theories turn out to be nonlocal [| Recently, using certain dynamical as- 
sumptions, one of us and P. Kurzepa succeeded to perform bosonization of QED in 2-1-1 
dimensions [Q. The resulting bosonic theory was found to be local and strongly inter- 
acting. As expected, in 3+1 dimensions the corresponding construction is significantly 
more complicated. In this letter we perform a first step of the analogous program in 3+1 
dimension: the construction of fermionic operators in terms of bosonic fields. 

We consider quantum electrodynamics with one Dirac fermion in 3+1 dimensions. 
Our basic philosophy and motivation to consider this particular theory are basically same 
as in 2+1 dimensions. They are described in detail in ref. |Q and need not be repeated 
here. Apart of obvious specific 3+1 dimensional features, like a nonabelian representation 
of the rotation group and different discrete symmetries, note that now the axial anomaly 
is present. This is different from 2+1 but reminiscent of 1+1 QED. We elaborate later on 
its implementation. We will work in the Hamiltonian formalism, and find the temporal 
gauge {Aq = 0) to be most convenient to our purposes. The Hamiltonian is 

H=^E^ + ifi' + H'i^d, + eAi)ij + m^^ (1) 

It is supplemented by Gauss's constraint, 

diEi = et/j^tP (2) 

We take the same approach as in ref.|Q. We would like to construct gauge invariant 
fermionic fields in terms of bosonic fields Ei and their conjugate momenta, which would 
solve the constraint eq. |^. In principle, as a next step we would like to calculate fermionic 
bilinears, including the components of the energy momentum tensor, in terms of these 
bosonic fields, and thereby bosonise the theory. In this paper we will only perform the 
first step of this program, namely construct the fermionic operators. 

The fields ipa that appear in eq.|l| are not gauge invariant. We therefore define 

^^^(x)=^„(x)<l>(x) (3) 

where, 

and ^ 

= (5) 

47r xr 



^The term "bosonization" is sometimes understood as obtaining an effective iow energy description 
of fermionic tlieories in terms of bosonic fields. This in many cases is a straightforward excersise [^. 
This is not what we mean by bosonization in this paper. Here we have in mind an exact mapping 
between the fermionic and bosonic Hilbert spaces, which should give an exact bosonic description of a 
fermionic theory at all energy scales. 
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is the electric field of a point like charge. The operator defined this way is gauge invariant, 
and coincides with the original fermionic field in the Coulomb gauge, in which the phase 
factor <l>(x) in eq.^ becomes unity. 

Our ansatz for the bosonic form of the fermionic operator is inspired by the results 
of ref. 0. In 2+1 dimensions the fermionic operators had the following basic structure]^: 

^2+i(x) = lim / deT{e)V^^\x + e)^x)V*^^\x - e) (6) 

The different factors in this expression have a clear physical meaning. The operator $ is 
the exact 2 + 1 dimensional analog of the operator defined in eq.^ It creates the electric 
field of a point electric charge at the point x. The role of this factor is to ensure the 
correct transformation properties of under the global electric charge transformation. 
The operator V{x) {V*{x)) is the operator that creates a magnetic vortex with the 
elementary fiux 27r/e (— 27r/e). This is a local scalar field in the framework of QED2+1 
p. Its appearance in eq.^ is natural, since bosonization should be akin to duality 
transformation, and V{x) is in fact the dual field in the theory 0. The fact that it is 
the square root of V that appears in the equation, makes the operator ip double valued. 
This is a necessary property, since otherwise ip{x) could not have the correct spinorial 
rotational properties. It is also in agreement with general arguments of ref. 0. V* also 
necessarily appears in the same expression, since should not carry magnetic fiux, and 
therefore the fiux that is created by V is cancelled by V*. This particular combination 
also insures anticommutativity of iIj{x). To make the product of local operators appearing 
in the definition of well defined, one has to point split V and $. The integration over 
the direction of the point splitting e with the appropriate phase factor r(e) projects the 
composite operator ip onto the spin 1/2 representation of the rotation group. 

To generalize this construction we should first understand what is the generalization 
of the vortex operator V to 3+1 dimensions. The operator V now should create an 
infinitesimally thin line of magnetic fiux. It will depend therefore on the curve. The 
explicit construction has been given in refs.|^, 0. 

V{C) = exp 1^ I d^y[a,{y - x,)E,{y) + b{x, - y)dMy)]] (7) 

Here aj(x) is a vector potential of an infinitesimally thin magnetic vortex along the curve 
C : eijkdjak{x) = 27rrj(x, C)6p^j.p{x), where rj(x, C) is the unit tangent vector to the curve 
C at point x, and 5perp(^) is the two dimensional delta function in the plane normal to Xj. 
The function b{x) satisfies di[b{x)]mod2TT = CLi{x). Since ai{x) has a nonvanishing curl, the 
function b{x) must have a surface of singularities ending at the curve C. For example, 
for a straight line C running along the X3 axis one has 

X ■ 

ai{x) = eij n ' n , i = l,2; 03(0;) = 

b{x) = e(x) (8) 

with the polar angle in the xi — X2 plane (Fig.l.). An alternative representation of V 

^In this formula we have not made explicit the difference between the first and the second components 
of the Dirac spinor. The details of the construction are given in ref. B. 
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Figure 1: The function Q{x). 
Figure 2: The standard set of curves C. 

is obtained by integrating by parts the diEi term in the exponentiah 

V{C) = exp{— / d^s'E, + -i <fs\b{x, - y)E,)} (9) 

where the integration is over the surface S of discontinuities of the function 6, which 
has the curve C as its boundary, and the second term in the exponential is the integral 
over the surface at spatial infinity. Note, that even though the surface S enters formally 
the definition of V , the operator actually does not depend on it, but depends rather 
only on the boundary C. This is the consequence of the fact, that changing the surface 
changes the phase by 27r/e / ds^Ei, which is an integer multiple of 2tt due to quantization 
of electric charge . 

Let us now define a standard set of operators Vr^{x, h), which are associated with the 
set of curves in Fig. 2. Those consist of half a circle with the center at the point x and 
radius t] in the plane specified by the unit normal vector n, and two semiinfinite segments 
of a straight line in the direction of the third axis. Note, that V*{x, —n) is obtained 
from V^(x, n) by a three dimensional parity transformation. Note, also that the product 
V^(x, h)V*{x, —n), is an operator that creates a unit strength closed t'Hooft loop, in the 
plane normal to n, centered at x with radius t]. 

We are now going to generalize the 2 + 1 dimensional construction of the Fermi opera- 
tors, by basically changing the product of a vortex and antivortex into an infinitesimally 
small t'Hooft loop of half integer strength. Consider therefore the following operators 

X{x)n = kA^^^e-^^ I '''y^^''^V^/\x,nMx)V;'/^{x, -n) (10) 

The first exponential factor here is the same as in 2 + 1 dimensions |Q. We have also 
introduced the ultraviolet cutoff A so that x has a dimension of the canonical fermi field. 
The constant k in eq. |l^ is should be chosen so as to normalize correctly the fermionic 
bilinears (see Q]). The operator $(x) is defined in eq. ^jin terms of an operator Ai. This 
operator Ai should be understood as a variable conjugate to Ei 

[E,{x),A,{y)]=z6\x-y)6,, (11) 

We also require, that eijkdjA^ = Bi. In these respects Ai{x) is similar to the original 
vector potential. One should remember, however, that since we are working on the gauge 
invariant subspace of QED4, Ai is a gauge invariant operator. In this sence, it can be 
thought of as a vector potential in a particular gauge. Up to now, we did not need to 
specify this gauge condition. We will return to this point later. 

■^We pick the third direction arbitrarily, for definiteness. As will become clear later, none of our 
calculations depend on it. 
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The operator y^/^ (y*^/^) obviously has different locahty properties than the operator 
V. Whereas Vc has support only on the curve C, and the dicontinuity surface of the 
function b{x) is unphysical, this is not quite true for V^^"^ (V*^^"^). The latter operator 
is double valued. By moving the surface of discontinuities from Si to 5*2, we multiply 
V^^"^ by a phase factor exp{^ !fsi-S2 '^^^^i}^ which due to quantization of electric charge 
can be either 1 or —1. Fixing the surface of discontinuities, is therefore equivalent to 
specifying what branch of the square root we choose in the definition of V^^"^. 

Although we do not have to specify this choice, and can just think about V^^'^ as a 
single valued operator defined on a double cover, it is sometimes convenient to visualize 
it as having a definite surface of discontinuities. We will always think about this surface 
as being parallel to the [1,3] plane. Note, that it is still true that V*^^'^{x, —h) the and 
V^/'^{x^n) are related by a parity transformation, and therefore the surface of disconti- 
nuities associated with the former operator, is also the parity transform of the surface 
of discontinuities associated with the latter. Therefore, the surface of discontinuitites 
associated with the closed loop, V^/'^{x^ n)V*^/'^{x, —n) in the limit of small rf is just the 
whole plane [1, 3]. This doublevaluedness, or effective nonlocality in the definition of V^^'^ 
is directly responsible for the anticommutativity of the operators defined in eq. Using 
the fact that for the function Q[x) defined in eq. || satisfies Q{x — y) — Q{y — x) = ±tt 
for any two points x and y, we find: 

Xnix)xn'iy)\x-y\»'n = Xn' iy)xhix) exp{i{e{x - y) - e{y - x))} = -Xrviy)Xnix) (12) 

The anticommutation relations between x{^) xKv) ^l^o follow immediatelly for 
X ^ y. In fact, the canonical anticommutator at coincident points is obtained, quite 
generally, only in the limit where the regulator in the definition of fermionic operators 
is removed t] ^ 0. This is also true in other dimensionalities [^]). In this limit the 
operator x does not depend on the direction fi and one obtains 

x{x)xHx) + xHx)x{x) = 2k^A^ (13) 

which for an appropriate choice of k becomes a 5 - function in the continuum limit. 

Complete analogy with 2 + 1 dimensions would require at this point to project the 
operators x onto appropriate representation of the rotation group. However, before doing 
that we have to address one more question, which was not present in 2 + 1 d. This is 
the action on the fermionic fields of the axial Ua{^) transformation. To implement the 
correct action of ^7a(1), we need to know the form of the axial charge in terms of bosonic 
variables. We infer this from the axial anomaly equation of QED^. 

= -^e.^x^F^'-F'^ (14) 



We take therefore the bosonic representation of the axial current as: 

4 = -^e'^'^'^^-^A. (15) 
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and in particular 

^5 = -^2^^B. (16) 

This is completetly analogous to 1 + 1 dimensional case. We want to stress, that eq. 
should not be thought of as a guess, but rather as the equation defining what we mean 
by the operator Ai. Up to now, we have only required that Ai is a canonical conjugate 
of the electric field Ei, and that eijkdjAk = Bi. Eq. |16| defines also the longitudinal part 
of Ai in terms of the gauge invariant quantity Jg , in a way consistent with the anomaly 
equation. This equation therefore can be understood as the gauge condition referred to 
earlier. 

Let us see what are the transformation properties of the operators x under the axial 
transformation. The operator $(a;), clearly commutes with Jg, since it only involves the 
momentum A. The relevant commutator is, therefore 

giaQYV2(a.^ h)v;^l\x, -n)e-*"«' = V^/\x, n)V;^'\x, -n)e*^"'^^(") (17) 

Here (j)n{x) is the magnetic fiux through the surface bounded by the circle of radius 
7] with the center in x lying in the plane with normal n. In other words it is mag- 
netic fiux through the area bounded by the closed t'Hooft loop created by the operator 
V^'^{x,n)V;^/^{x,-n). 

We can now define eigenoperators of the axial charge with eigenvalues ±1 

Xn(x)+ = j dae'''^' Xn{x)e-"'^" e-"" = J daXn(a;)e*"(^^'^(^)-^) 

Xn{xy = J d«e*"^'xn(x)e-*"^'e'° = J daXn{x)e"'^^^^^'''^+^^ (18) 

The integral over the transformation parameter a can be now easily performed, and 
we find 

Xh{.xy = Xhix)6{(f)n{x) - 1) 

Xn{x)- = Xn{x)6{Mx) + I) (19) 

This result has a very natural interpretation. The operator Xn{x)^ is nontrivial only 
on states which contain a unit magnetic fiux linked with the t'Hooft loop in the definition 
of Xh{x)~^- Note, that the extra factor does not change the anticommutation relations 
of the operators. This is obvious, since this extra piece measures magnetic field at the 
location of the operator xix), whereas the operator xiv) creates only a closed loop of 
fiux at its location. More precisely 

e'''^^^-^^''^V^/\y, n')V;^/\y, -n') = V^/\y, n')V;^/\y, -n')e*"^^^(")e*t'(f^i'f^2) (20) 

where l{Ci, C2) is the linking number of two small circles with radius 77 centered at points 
X and y with normals n and n' repsectively. 
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We now project these operators onto spin 1/2 representations of the rotation group. 
We use the chiral representation of the Dirac 7 matrices, in which each chirahty com- 
ponent transforms under rotations independently. The complete Dirac multiplet is then 
constructed as follows: 



i)i{x) = / dnYi/2{h)xhixy 
i^six) = J dnYi/2in)xn{xy 

^p^{x) = J dnY^i/2{n)xhix)~ (21) 
Here ^±1/2(5^) are monopole harmonics 

Y±y2{<P,0) = Te^HlTcos9Y/' (22) 

By construction therefore ipi,ip2 constitute a rotational doublet of chirality +1, and 
ipSjipA - a rotational doublet of chirality —1. 

Eq.^ is our final result. To conclude, we have constructed in the framework of 
QED4 the four dimensional multiplet of fermionic operators entirely in terms of bosonic 
fields Ei and their conjugate momenta Ai. The fermionic operators have canonical an- 
ticommutation relations in the limit of zero regulator, carry unit electric charge, can 
be decompose into two representations of the axial charge with axial charges ±1 and 
have correct transformation properties under rotations. Note also, that under the parity 
transformation 0(n) —(j){n), and Xn X-h- Under the parity transformation there- 
fore the rotational multiplets ipi,ip2 and ip^iil)^ are transformed into each other, which 
is precisely the right parity transformation properties of Dirac fermions in the so called 
chiral basis of 7 matrices |Tl|p|. 



The bosonic form of the fermionic operators has a very clear intuitive interpretation. 
The operator creates the electric field of a pointlike electric charge and the magnetic 
field of an infinitesimally small half integer magnetic fiuxon. The extra factor required 
to represent the axial charge, vanishes unless the operator acts on a state which has a 
unit magnetic fiuxon with a linking number ±1 relative to the fiuxon created by the 
fermionic field itself. The next logical step is to calculate fermionic bilinears in terms 
of the bosonic fields. This calculation was rather involved already in 2 + 1 dimensions, 
and is even more complicated in the present case. Work along these lines is currently in 
progress 
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*We wish to make a remark on the parity transformation properties of the operators. If we were to 
define the action of parity on y^/^ as V'^/'^{n) V*'^/'^{-n) it would follow that x ^ e*"^ / ''"^'^'^^^x 
which would not give the correct parity transformation properties. However, e'~ -f ^'■^'(y) jg a,n operator 
with eigenvalues ±1. Since x is a double valued operator, we are free to change our definition of parity 
so that X — >• X- This does not change the parity transformation properties of all single valued operators 
and is therefore perfectly consistent with the standard parity transformations of Ai and Ei. 
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